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g : ABSTRACT 

(N 

PLi, Active Galactic Nuclei can produce extremely powerful jets. While tightly 

"*^ I collimated, the scale of these jets and the stellar density at galactic centers implies 

OO I that there will be many jet /star interactions, which can mass-load the jet through 

stellar winds. Previous work employed modest wind mass outflow rates, but this 

^.^ _ does not apply when mass loading is provided by a small number of high mass- 

lO ! loss stars. We construct a framework for jet mass-loading by stellar winds for a 



t 



broader spectrum of wind mass-loss rates than has been previously considered. 



OO 

in 

"^ I Given the observed stellar mass distributions in galactic centers, we find that 

\^ [ even highly efficient (0.1 Eddington luminosity) jets from supermassive black 

^ I holes of masses Mbh ^ W^Mq are rapidly mass loaded and quenched by stellar 

i-^ ' winds. For IO^Mq < Mbh < lO^M©, the quenching length of highly efficient jets 

I I is independent of the jet's mechanical luminosity. Stellar wind mass-loading is 

H I unable to quench efficient jets from more massive engines, but can account for the 

52 I observed truncation of the inefficient M87 jet, and implies a baryon dominated 

composition on scales > 2 kpc therein even if the jet is initially pair plasma 

dominated. 



Subject headings: stars: mass-loss-stars: winds, outfiows-lSM: bubbles-ISM: jets 
and outflows-galaxies: active-galaxies: jets 



1. Introduction 

Active Galactic Nuclei (AGN) can produce extremely energetic jets. While these jets 
can be very tightly collimated out to kilo-parsec scales, the volume contained by these jets is 
large and they will inevitably encounter many stars while passing through the dense centers 
of galaxies. Komissarov (1994) established stellar winds as a source of jet mass-loading. 
He modeled tails blown downstream but neglected the finite size of the initial wind-blown 
bubble. Komissarov found that for a total mass loading rate supplied by individual stars all 
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supplying mass at the averaged mass-loss rate of IO'^'^Mq/jt, the tails merged and spread 
their mass throughout the entire jet cross section. He also found that when the same total 
mass-loading rate with a fixed average mass loading rate is provided by a small fraction of 
high mass-loss stars in the jet, their tails would neither merge nor cover the jet cross section. 

Here however, we find that bubbles blown by stellar winds inside a jet can quench the 
jet irrespective of the tails. AGB stars and massive main sequence stars can have mass-loss 
rates as high as 10~^MQ/yr, and the bubbles blown by even a small number of such stars can 
block substantial portions of the jet (Bednarek & Protheroe 1997). We examine the effect 
of bubbles blown by a variety of stellar winds and show that while high mass, low terminal 
velocity winds (such as those from AGB stars) are poorly spread throughout the jet, such 
winds are nevertheless strong enough to significantly affect jet dynamics. Bubbles efficiently 
spread the material of winds launched by stars with low mass-loss rates or high terminal 
velocities throughout the jet. 

In Sec. 2 we derive the conditions that determine whether stellar winds can truncate 
a given AGN jet. These depend on the total stellar wind mass-loss rate encountered by 
the jet and the fraction of the jet cross section intercepted by each wind. In Sec. 3 we 
discuss parameter regimes of jet mass loading for application to specific cases. In Sec. 4 we 
apply the formalism to the kilo-parsec scale jet of M87, and find that our estimates for the 
mass-loading rates required to stop the jet agree with the total amount of stellar mass-loss 
contained within the jet. We conclude in Sec. 5. 



2. Conditions for jet mass truncation via stellar winds 

To study the stellar wind mass-loading of a relativistic jet we must compare the total 
mass-loading rate Mt of all wind sources to the jet's mechanical luminosity Lj (hereafter, 
reference to "jet luminosity" will imply the total mechanical luminosity). Taking 7 as the 
jet's Lorentz factor, we see that if '-jMnp(? < Lj, the mass-loading is too weak to affect the 
jet while if 'jM^c^ > Lj the jet can be slowed or stopped. We define 

Mj = Ljhc^ (1) 

as the critical mass-loading rate. If a black hole of mass Mbh and corresponding Eddington 
luminosity Lem launches a jet of mechanical luminosity Q.ILemi Mj ~ {2/'~^){Mbh/Mq) x 
IQ-^^Mq/ji. For Mbh = W^Mq and 7 = 2 then, Mj = W-^Mq/jt. This value of Mj is 
approximately the upper limit for known stellar winds, implying that such a jet could be 
impeded by a single object. However, the condition Mt = Mj does not suffice to slow the 
jet. A second requirement is that the mass be spread throughout the jet. We express this 
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additional condition by considering the fractions / of the jet's cross-sectional area intercepted 
by stellar wind bubbles. 

The two conditions required for quenching the jet are then 

^ M^ = Mt>Mj, (2) 

stars n 

and 

E fr.>h (3) 

stars n 

where M„ and /„ are the wind mass loss rate and wind covering fraction for the nth star. 
In Eq. (3), we have ignored the possibility of two bubbles intercepting the same arc of the 
jet and the shielding of one stellar wind by another. 

The fraction of the jet blocked by a wind will approach a maximum J'm = f{t = oo) 
that we will show to be j'm = M/Mj in Sec. 2.1. This is precisely the fraction of the jet 
that the wind can load according to Eq. 2, and thus winds characterised by / = /m may 
be considered to have maximally loaded the jet. We define an interception function for each 

star 

/ = f/fM (4) 

and use it to rewrite Eq. 3 as 

^"^ r \^ T f l^stars n -'"^^" ^ i /|-N 

In — / ^ ^njM,n — "v — -•-• I'^j 



stars n stars n 



Mj 



Since / < 1, any system meeting the mass- mixing requirement (Eq. 5) also meets the total 
mass-loss requirement (Eq. 2) and hence Eq. (5) becomes the necessary and sufficient 
requirement for jet truncation via stellar winds. 

Eq. (5) allows us to relate the volume of a mass-loaded jet (through the number of stars 
contained therein and the stellar number density), the type of the stars (through / and M) 
and the luminosity of the jet (through / and Mj). We explicitly derive I in the next section 
and outline an estimate for Mt in section 2.2. In section 3.3, we will see that I depends 
linearly on Mj for some star /jet combinations. We can see that if such stars dominate jet 
mass-loading, the jet's luminosity drops out of Eq. (5). In these circumstances, the term to 
the left of the inequality in Eq. (5) depends only on Mt-, as will be discussed in section 3.4. 



2.1. The interception function 

As shown schematically in Fig. 1, the bubble blown by a wind has radius r and intercepts 
a fraction / = r'^ /R^j of the jet. We characterise a stellar wind by its mass-loss rate M and 
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terminal velocity v^. The star crosses the jet perpendicularly with a velocity Vc- If the jet 
is tightly collimated, then its radius Rj will be much less than its length. We presume that 
the jet's luminosity Lj is evenly distributed over its cross-sectional area Aj = vri?j. 

As long as the stellar winds are reasonably far from the central engine, a relativistic 
jet velocity will dominate the other velocities (thermal speeds, foo and Vc)- We consider a 
wind-blown bubble to be an onion-like set of concentric shells, of which only the outermost 
will interact with the jet. Once a shell of mass Am has absorbed energy '-fAmc^ from the 
jet, the shell will be carried away by the jets and thus lost from the bubble. This allows us 
to ignore lateral pressure, as any portion of the stellar wind affected by such forces would 
have already been carried downstream. Thus, the bubble will expand at its (unchanged) 
terminal velocity v^o while simultaneously losing its outer shells to the jet. In a time At, 
the wind absorbs an energy AE = Lj{r'^ / R^j)At from the jet and accordingly loses a shell 
of mass Am = AE/{'yc^) = Mj{r^/Rj)At. This shell has width Ar = {Am/M)v^ = 
{r"^/ fMR'j)vooAt, where we have used /m = M/Mj (to be justified below). We can write 
the equation for the growth of the wind-bubble as follows: 



dr Ar / r^ 

'dt^''°°~Ai^'''^(7^R 



^oo - ^r: = ^oo ( 1 - J „2 ) • (6) 



The above equation can be solved to give 

r(t) = Vf^Rj^;j^^y^ (7) 

where 

r]{t) = 2va/y^Rj (8) 

and 7f{t)/4: is the ratio of the jet mass flux the bubble would intercept, were it unaffected 
by the jet, to the stellar wind's mass loss flux. From Eq. (7) we find that 

As in the above equation /(oo) = j'm, our definition of /m is justified, and the interception 
function / of (4) becomes 
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2.2. The total mass loading from stellar winds and the role of the IMF 

In the absence knowing the specific stars in a jet, we can use an initial mass function 
N{m). In this case Eq. (2) becomes 

/rrif 
N{m)M{m)dm > Mj (11) 

szars n ^» 

and Eq. (5) becomes 

N{m)I{m)M{m)dm > 1. (12) 

star*- 7t -^« 

In general, I{m) will depend on the age of the stellar population in question. In Sec. 3 we 
find regimes where / can be approximated. Eq. (11) is simpler as it does not invoke I. 

In estimating the relevant mass loss M^ of the stars contained in a jet, we first note that 
although SN deposit a lot of energy, they are not a significant contributor to stellar mass- 
loss for typical Galactic supernovae rates of 1-2 per century. We therefore assume that the 
averaged total mass-loss rate of a population of stars closely approximates its averaged wind 
mass-loss rate. We also assume that a star loses a fraction a ~ 1 of its mass over its lifetime 
r. The lifetime of a star behaves as r ~ em/ L where m and L are its mass and luminosity 
and e = TqLq/Mq provides appropriate scaling. Accordingly, over their entire lifetimes, 
stars exhibit average mass-loss rates of M ~ aM/r = aL/e. While very low mass-stars lose 
only a small fraction of their mass over galactic lifetimes, their aggregate luminosities and 
mass-losses are negligible. 

For an arbitrary IMF N{m), the total stellar luminosity and total averaged mass-loss 
rates from a population of stars are then 

L{m)N{m)dm (13) 



and 

Mt= I a—N{m)dm = / -L{m)N{m)dm = -Lt = a^— ^ a x 10"^°— Me/yr. 

Jnii T Jmi £ ^ Tq Lq Lq 

(14) 
Remarkably, N{m) does not affect the final form of the above equation and Mt depends only 
on Lrp- Applying Eq. (14) to the Milky Way, we find a stellar mass- loss rate Mt ~ 4M0/yr, 
which is comparable to estimates for the observed star formation rate (e.g. Diehl et al. 
(2006)). This justifies use of Eq. (14) for the total stellar mass-loss rate. 
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Despite not affecting Eq. (14), the form of A^(m) can influence tlie effect that stellar 
wind mass-loading has on a jet through / (via Voo and M). In addition, the total mass loss 
can be dominated by a small number of massive winds. If this number is small enough, the 
process of estimating M^ by averaging and summing over the winds contained in the jet over 
long periods of time will not accurately describe Mt at any given instant. If the condition 
given by Eq. (5) is on average marginally fulfilled for the characteristic jet length, significant 
fluctuations of the stellar population summed over in Eq. (5) can create significant jet length 
fluctuations. 



3. Parameter Regimes for Analysing Wind Mass Loading and Jet Truncation 

While Eq. (5) provides the fundamental condition for whether stellar wind mass-loading 
can truncate jets, its implementation is complicated by the need to accurately assess the wind 
parameters of the constituent stars and individual jet/wind interactions. Here we examine 
some specific ranges that parameters in Eq. (5) can take in order to understand the effects 
of individual jet/wind interactions as well as to find regimes in which global effects can be 
easily calculated. 



3.1. Strong winds (/m > 1) 

For a jet containing a stellar wind with Jm > 1, we can use Eq. (7) to calculate the time 
tfe it will take for the associated wind bubble to grow and intercept the entire jet. Setting 
r = Rj in Eq. (7) gives 

t,^^MlLn(l±y^\. (15) 

As /m increases in the regime /m > ^,tb rapidly approaches -Rj/foo, and the bubble expands 
essentially unimpeded. If f oo > Vc, the star will eventually intercept the entire jet. If instead 
■^oo < Vc the wind will have insufficient time to intercept the entire jet and must be treated 
as winds with Jm < 1, which we discuss next. 
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3.2. Weak winds (/m < 1) 

Here we examine tlie bubble after a time t = tc = Rj/vc required for the star to cross 
halfway through the jet. Accordingly, Eq. (8) becomes 

Vc^v{tc) = ^ — . (16) 

Substituting this into Eq. (10) gives the characteristic Ic for weak winds: 

(17) 



ei- + 1 

In Fig. 2 we plot Ic as a function of ?7c- If Vc is large, / approaches its maximum value /m- 
This will occur for weak (/m <^ 1) and fast (voo ^ fc) winds. 



3.3. High Mixing vs. Low Mixing Winds 

There is a wide range of possible wind and jet parameters (e.g. M can vary by 10 orders 
of magnitude) but useful conclusions can be drawn merely by dividing the wind-jet systems 
into "high mixing" [rjc > 1) and "low mixing" {rjc < 1) categories. For "high mixing" winds, 
/ ~ 1 and Eq. (5) simplifies to Eq. (2). This "high mixing" category includes both low 
mass winds (/m <^ 1) and high mass-loss, fast winds (foo ^ Vc) such as Wolf-Rayet stars, 
which can single-handedly slow an entire jet. 

Low mixing winds {rjc < 1), such as those from AGB stars, have high mass loss rates 
but low terminal velocities. As seen in Fig. 2, for such winds 

/(„)-f = f f^)'. (18) 

4 Im \Vc J 

Unless the wind-jet interaction takes place close to a galactic center, the galaxy's velocity 

dispersion a will approximate the star's crossing velocity Vf.. Accordingly we rewrite Eq. 

(18) as 

1 



/c=^m . (19) 

An interesting implication of Eq. (19) is that the fraction of the jet blocked by a low mixing 
wind, / = r'^/R'j = IcI'm = {voo/o')'^, is independent of both the wind's mass-loss rate and 
the jet's mechanical luminosity. This is because the wind-blown bubble stays small enough 
that the second term on the right of Eq. (6) remains negligible. 



Because M = 10~^Mq/jt is a reasonable upper limit on wind strengths for both fast and 
slow wind stars, Eq. (16) demonstrates that all stars are "high mixing" for any adequately 
powerful {Mj > lO~^M0/yr) jet. However, stellar winds do not provide adequate mass- 
loading rates to substantially slow such jets. Stellar winds are able to slow weaker jets, but 
the broad range of values that I can take for "low mixing" winds that enter the jet requires 
careful modeling to use Eq. (5) in a meaningful way. 



3.4. Jet luminosity regimes 

As seen above, the jet strength Mj does not influence Eq. (5) for "low mixing" winds. 
If we presume that ^ In^n is dominated by "low mixing" winds from AGB stars, Eq. (5)'s 
resulting lack of dependancy on Mj allows us to draw general conclusions from Eq. (5) 
without knowing the behaviour of / in detail. We use Eq. (18) to rewrite Eq. (5) as 

J I cliar_ \ J_ ^ ^ ^20) 

/ nMr ~ ^ 



M,har ^ Cr / gMj 

where Mchar and Vchar are the mass-loss rate and terminal velocity of the characteristic peak 
mass-loading winds. The parameter g approximates the inverse fraction of the total mass-loss 
contributed by such winds. Galactic stellar velocity dispersions satisfy a ~ 200 — 300km/s. 
For winds from AGB stars, Mchar ~ 10^^ — lO~^M0/yr and Vchar ~ lOkm/s. We therefore 
find that 

MT-6xlO-'g( ^Vf^^^ I ^'^""'', iMe/yr (21) 

^ V250km/sy VlOkm/sy \^lO-5M0/yry ^'^ ^ ' 

approximates the total mass-loss required to truncate such a jet. This calculation applies 
to jets with Mj < lO'^Mchar- For significantly weaker jets, the assumption that UlnMn is 
dominated by a narrow range of stars breaks down and the jet can be more easily quenched 
(although the quenching depends sensitively on /). Nevertheless, for an intermediate range 
of jet luminosities, Eq. (21) does apply. In this regime, the stellar mass-loss required to 
truncate the jets is independent of the jet's luminosity for mass-loading dominated by "low 
mixing" winds. 

We are now able to define three jet luminosity regimes, characterised by Mj. (1) Strong 
jets, with Mj > lO^^MQ/yr, experience all winds as "high mixing." In this case Eq. (5) 
can be approximated by Eq. (2). However, these jets are so strong that stellar winds do 
not represent a viable mechanism for their truncation. (2) Intermediate strength jets, with 
lO^^Mo/yr < Mj < lO^^M0/yr, may have their mass-loading dominated by "low mixing" 
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stars. If so, such jets will propagate until the total wind mass-loading rate into the jet is 
Mt ~ IQ-^Me/yr (Eq. 21). (3) Low luminosity jets with Mj < 10'^ Mq/jt are too weak 
to withstand stellar wind mass-loading for long, but the details of the loading depend more 
strongly on the stellar population than the other two categories. 



4. Application to observations 

As long as the total mass-loss of stars beyond the jet's truncation length remains signif- 
icant, stellar wind mass-loading is a good candidate for jet truncation. The stars enter the 
jet from the side and would likely induce an more extended but staccato truncation of the 
jet compared to jet truncation via ISM pile up. Both forms of mass loading would imply a 
baryonic jet composition at the truncation scale. Once mass loading reduces the longitudinal 
velocity of the jet to its internal sound speed, the jet will expand into a quasi-spherical lobe. 
In this respect, radio lobes are a natural prediction of jet truncation via mass loading. 

If jet truncation via stellar winds is assumed, then Eq. (5) can be used to constrain the 
jet properties and the stellar population it contains. In particular, by assuming the left side 
of Eq. (5) exceeds unity, we can use any two of Mj, Mt and / to calculate the third. 

Alternatively, one can use Eq. (5) to evaluate the effect that stellar wind mass-loading 
has on the jet by knowing the jet's total mechanical luminosity Mj, total mass-loading rate 
Mt and the contained stellar population profile (which sets I). In this case, the inequality 
is assessed, not assumed. 



4.1. Application to M87 

The elliptical galaxy M87 has a well studied jet with luminosity Lj estimated to be 
between 10^^ and 10^^ erg/s (Biretta et al. (1991), Reynolds et al. (1996)) and a 7 of about 6 
(Biretta et al. (1999)). It follows that Mj > 3 x IQ-^Mq/jt. This value for Mj is interesting 
as it can be attained and indeed surpassed by single stars and suggests that the jet is in 
our intermediate luminosity regime (see Sec. 3.4). The collimated jet has a length of 2 kpc 
(Perlman et al. (2001)), and is therefore truncated well within the confines of M87. Stellar 
wind mass-loading is therefore a plausible cause of the jet's truncation (implying the jet 
could be longer lived that its length alone would indicate). 

Faber et al. (1997) provide values for the luminosity density near the center of M87. 
The jet contains a total stellar luminosity of roughly 2 x lO"^ Lq in the V-band and so, by Eq. 
(14), a stellar mass loss of Mt ~ 10^'^Mq/jt. This fits our estimate of the total mass-loss 



-10- 



rate Mt required to quench intermediate luminosity jets from Eq. (21), suggesting that "low 
mixing" stars dominate the mass-loading process. 

Reynolds et al. (1996) find evidence that the jet may be initially pair plasma dominated 
near the launch point, but this does not contradict a subsequent increase in baryon frac- 
tion due to mass-loading via stars on kpc scales. Reynolds et al. (1996) also estimate the 
mechanical luminosity of the M87 jet via its infiuence on the radio lobes. This estimate is 
independent of the jet composition. 



5. Conclusion 

We have modelled jet mass-loading via stellar winds by focusing on the bubbles blown 
by these winds. We have derived a simple condition between the jet and wind parameters 
that determines whether mass-loading and can truncate jet propagation. From this relation, 
jet truncation via stellar winds can be used to constrain the jet's mechanical luminosity 
(and baryon fraction), the total mass-loss rate of the stars contained within the jet, or the 
population profile of those stars. 

We have applied our model to the kilo-parsec scale jet of M87, and find that jet trun- 
cation via stellar winds provides a consistent interpretation of the length of that jet. This 
would also imply that stellar winds at the center of that galaxy are dominated by high 
mass, low velocity winds, and that the jet is baryon dominated at the scale of its truncation, 
regardless of its initial launch composition. 
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Energetics. EGB acknowledges support from NSF grants AST-0406799, AST-0406823, and 
NASA grant ATP04-0000-0016 (NNG05GH61G). 
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Fig. 1. — Schematic representation of a star and wind in our model. The star is crossing 
the jet, of radius Rj, with velocity Vc- The wind, of radius r is expanding at its terminal 
velocity v^o- 
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Fig. 2. — Interception function (/ = f / fu) as a function of r^c = 2t>oo/ ("WcVTm) (solid curve) 
as well as its low rjc limit rJl/A (dashed curve). A value of 1 implies that the bubble blown 
by the stellar wind would, having crossed halfway through the jet, have expanded to the size 
at which it intercepts energy flux from the jet equal to that produced by the star's mass-loss 



